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ABSTRACT 

We introduce new methods for robust high-precision photometry from well-sampled images 
of a non-crowded field with a strongly varying point-spread function. For this work, we used 
archival imaging data of the open cluster M37 taken by MMT 6.5m telescope. We find that 
the archival light curves from the original image subtraction procedure exhibit many unusual 
outliers, and more than 20% of data get rejected by the simple filtering algorithm adopted by 
early analysis. In order to achieve better photometric precisions and also to utilize all available 
data, the entire imaging database was re-analyzed with our time-series photometry technique 
(Multi-aperture Indexing Photometry) and a set of sophisticated calibration procedures. The 
merit of this approach is as follows: we find an optimal aperture for each star with a maximum 
signal-to-noise ratio, and also treat peculiar situations where photometry returns misleading 
information with more optimal photometric index. We also adopt photometric de-trending based 
on a hierarchical clustering method, which is a very useful tool in removing systematics from 
light curves. Our method removes systematic variations that are shared by light curves of nearby 
stars, while true variabilities are preserved. Consequently, our method utilizes nearly 100% of 
available data and reduce the rms scatter several times smaller than archival light curves for 
brighter stars. This new data set gives a rare opportunity to explore different types of variability 
of short (^minutes) and long (^1 month) time scales in open cluster stars. 

Subject headings: methods: data analysis — open clusters and associations: individual (M37) — stars: 
variables: general — techniques: photometric 


1. INTRODUCTION 

We are poised on the threshold of unprece¬ 
dented technical growth in wide-held time domain 
astronomy, where ground-based observations yield 
very precise measurements of stellar brightness 
from high-volume data streams. So far, wide- 
held time-series surveys has been spearheaded 
by relatively small telescopes since they are sup- 
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ported by large held of view (FOV) instruments 
opera ting with high duty cycle (see iBecker et al.l 


20041 for a summary of optical variability surveys). 
Within the last decade, the advent of large mo¬ 
saic CCDs has facilitated the coverage of large sky 
area even f or large-aperture tele scopes (e.g., MMT 
Megacam: McLeod et ^ _ 2000t ESQ Very Large 


Telescope Omegacam: Kuiiken et al.l l2002t Sub¬ 


aru Suprime-Cam: Miyazaki et al.l 200^ CHFT 
Mega cam: iBoulade et al.l 20031 : iPTF: Kulkarnil 
2OI3II . Although these facilities are generally 
devoted to imaging surveys, researchers are at¬ 
tempting to utilize them for short- and long-term 
variab ility surveys wi t h shor t -cadence exposures 

20091: 


Hartman et al. 2008ai iPietrukowicz et al 


Randal et al .11 20 111) . Such wide-held imag¬ 


ing systems haye enabled us to obserye hundred 


1 










































of thousands of target stars simultaneously and 
also to detect various variability phenomena. A 
remarkable thing about these surveys is that the 
fraction of variable sources increases as the photo¬ 
metric precision of the survey improves. For this 
reason, it is important to improve the accuracy in 
photometry. 

Another key issue in wide-held time-series pho¬ 
tometry is the removal of temporal systematics 
from a single image frame or several consecutive 
image frames. It has recently become known that 
systematic trends in time-series data can be dif¬ 
ferent and localized within the image frame when 
the FOV is large. Such spatially localized patterns 
may be related to subtle point spread function 
(PSF) differences and sky condition within the de- 


tector FOV (e.g., Kovacs et aLll2005: 

^enner et al. 

I 2 OO 8 I: iBianco et al. 20091: Kim et al. 

200911. As 


these patterns change in time, we can see how 
the temporal variations of systematic trends af¬ 
fect the brightness and shape of light curves di¬ 
rectly. The time-scale of systematic variation is 
sometimes comparable to short-term variability, 
such as transits or eclipses, and in some cases even 
long-term variability. Thus, it is often difficult to 
identify and characterize true variabilities. 

In this paper, we introduce a new photometry 
procedure, called multi-aperture indexing, which 
is suited to analyzing well-sampled wide-held im¬ 
ages of non-crowded helds with a highly vary¬ 
ing PSF, such as those produced by wide-held 
mosaic imagers on large telescopes. We apply 
this procedure to archival imaging data from the 
MMT/Megacam transit survey of the open clus¬ 
ter M37 (Hartman et al. 2008a), demonstrating 
a substantial improvement over the existing pho¬ 
tometry. Section 2 describes the MMT imaging 
database and identihes problems in the existing 
photometry which motivated the development of 
our new methods. Section 3 describes the multi¬ 
aperture photometry that utilize newly dehned 
contamination index and carefully tuned calibra¬ 
tion procedures, including the results of the basic 
tests to validate our approach. Section 4 gives 
an in-depth discussion about systematic trends 
in time-series data and suggests an efficient way 
for identifying, measuring, and removing spatio- 
temporal trends. Section 5 describes the effects of 
new calibration on period search, and we summa¬ 
rize our main results in the last section. 


2. ARCHIVAL MMT TIME-SERIES DATA 
OF M37 


Hartman et al.l (|2008all have conducted a study 


to hnd Neptune-sized planets transiting solar-like 
stars in the rich open cluster M37. The observ¬ 
ing strategy was carefully designed for a transiting 
planet search by several considerations (e.g., the 
reliability of exposure time per frame, the effects of 
pixel-to-pixel sensitivity variations, and sensitivity 
of filter). Their work did not reveal any transiting 
planets, but it did provide a rare opportunity to 
explore photometric variability at relatively high 


tempor al resolution with 30-90 s. iHartman et al 


( 2008bll discovered 1430 new variable stars, includ¬ 
ing very short-period eclipsing binaries (e.g., V37, 
V706, VI160) and 6 Set-type pulsating stars (e.g., 
V397, V744, V1412). 

We used the same data set on the open clus¬ 
ter M37. A detailed discussion of the observa¬ 
tions, original data red uction, and light curve p ro- 
duction is described in Hartman et al.l (j2008alft3 b 
The data archive consists of approximately 5000 
r'-filter images taken over 24 nights with the wide- 
held mosaic imager (Megacam) mounted at the 
//5 Cassegrain focus of the 6.5m MMT telescope. 
Note that Megacam is made up of 36 2048 x 4608 
pixel CCD chi ps in a 9 x 4 patte rn, covering a 
24'x24' FOV (iMcLeod et al.ll2000ll . This instru¬ 
ment has an unbinned pixel scale of 0'b08, but it 
was used in 2 x 2 binning mode for readout. 

The observation l ogs are summarized in Table 
1 of Hartman et al. ( 2008all . In brief, the r'-band 
time-series observations were undertaken between 
2005 December 21 and 2006 January 21, with a 
median FWHM of 0.89 ± 0.39 arcsec. Exposure 
times are chosen to keep an r ~ 15 mag star as 
close to the saturation limit, which is expressed as 
a function of seeing conditions. With an average 
seeing ~0".89 on images, the quality of the images 
is good to achieve high-precision light curves (less 
than 1% rms value) down to 20. In addition to 
the imaging data set, this database includes light 
curve data sets for a total of 23,790 sources de¬ 
tected in a co-added reference image. Theses light 
curves are obtained by the image subtraction tech- 
nique using a modified version of ISIS software 
( Alard fc Lupton 1998t lAlard 2000ll . 

As shown in Figure [U however, the raw light 
curves from the original image subtraction pro- 
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Fig. 1.— Example of MMT light curves for the 
brighter stars in archive. The dashed lines are 
weighted spline approximation with ±3 -(t control 
limits (dotted lines). These light curves contain 
outlier points that significantly increase the rms 
scatter of the raw light curves. 


cedures exhibit many unusual outliers, and more 
than ~15% of data get rejected by a simple filter¬ 
ing algorithm after cleaning procedures. In prac¬ 
tice, brutal filtering that is often applied to remove 
outlying data points can result in the loss of vi¬ 
tal data, with seriously negative impact to short¬ 
term variations such as flares and deep eclipses. 
We also find that the image subtraction technique 
often resulted in measurement failures from sev¬ 
eral frames due to poor seeing or tracking prob¬ 
lem. After removing these bad frames, it leads to 
loss of additional ~5% data points from most light 
curves. In order to overcome this problem, we have 
re-processed the entire image database with new 
photometric reduction procedures. 


3. NEW PHOTOMETRIC REDUCTION 


3.1. Preparation for Photometry 


We followed the standard CCD reduction pro¬ 
cedures of the bias correction, overscan trimming, 
dark correction , and fl at-fielding as described in 
Hartman et al. (2008a). The individual CCD 
frames were calibrated in IRAF, using the mo¬ 
saic data reduction package MEGARED0 The 


^The 64-bit version of Megacam re¬ 
duction package is available from 


first step is to correct the pixel-to-pixel zero-point 
differences that are usually described by the sum 
of a mean bias level and a bias structure. As 
the bias frames were not separately taken during 
the time of the observations, the mean bias level 
was subtracted from each image extension using 
an overscan correction and so we cannot remove 
any remaining bias structure from all overscan- 
subtracted data frames. According to descrip¬ 
tion in Matt Ashby’s Megacam reduction guideH 
the bias structure can be very significant in some 
small regions such as the portions of the arrays 
close to the readout leads. The dark currents are 
normally insignificant for Megacam so that correc¬ 
tions are not needed even for long exposures. The 
next step is to correct pixel-to-pixel variations in 
the sensitivity of the CCD; we used the program 
domegacamf lat^. This program determines the 
scaling factors to correct the gain difference be¬ 
tween the two amplifiers of each chip by finding 
the mode in the quotient of pixels to the left and 
right of the amplifier boundary, and then flattens 
each of the frames with a master flat field frame. 
It is worth mentioning that the sky conditions 
were rarely photometric during the observing run, 
with persistent light cirrus for most of the nights. 
Therefore it was only possible to obtain twilight 
sky flats on a handful of nights (dome flats were 
not possible). 

We removed bad pixels using the Megacam 
bad pixel masks distributed with the MEGARED 
package. The values of bad pixels are replaced 
with interpolated value of the surrounding pixels 
using the IRAF task f ixpix. The numerous single 
pixel events (cosmic-rays) were iden tified and re¬ 
move d using the LACosmic package ( van Dokkuml 
200 lb . 


The Megacam data already have a rough World 
Coordinates System (WCS) solution that is based 
on a single value of the telescope pointing. To 
update these with a more precise solution, we ap¬ 
plied astrometric correction to each CCD in the 


https: //www. cfa.harvard.edu/~bincleod/Megared/. 

^The detailed reduction procedures of MMT/Megacam 
are described by Matthew L. N. Ashby at 
https: //www. cf a. harvard. edu/~mashby/inegacain/inegacain_f rames .html 
^This C routine is written by J. D. Hartman to reduce I/O 
overhead. It reads in a list of Megacam mosaic images and 
a Megacam mosaic fiat-field image, and then writes out the 
flattened image over the existing image. 
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mosa ic using the WCSTools imwcs program (|Mink 


2002h . The new solution is derived by minimizing 
the differences between the R.A. and decl. posi¬ 
tions of sources in a single CCD chip and their po- 
sitions listed in the 2M ASS Point Source Catalog 
( Skrutskie et al.lBoOGll . The resulting astrometric 
accuracy is typically better than 0".l rms in both 
R.A. and decl. 


3.1.1. Building a Master Source Catalog 


Typically, a point or extended source detection 
algorithm is applied to each frame independently 
and it always requires criteria for what should 
be regarded as a true detection. In obtaining 
the pixel coordinates for all objects in the M37 
fields, this procedure often misses some objects 
when the detection threshold approaches the noise 
level. Also it needs a substantial effort to match 
the objects that are detected in only some of the 
frames. Our approach is as follows: a complete list 
of all objects is obtained from a co-added reference 
frame, and then the photometry is performed for 
each frame using the fixed positions of the sources 
detected on the reference. Since the relative cen¬ 
troid positions of all objects are the same for all 
frames in the time series, we can easily place an 
aperture on each target and measure the flux even 
for the stars at the faint magnitude end. 


We constructed the reference frame for each 
chip from the best seeing frames using the SWarjjl 
software. Benefiting from a highly accurate as¬ 
trometric calibration of input frames, we were 
able to improve the quality of co-added images. 
In the SWarp implementation, the pixels of each 
frame were resampled using the LanczosS convo¬ 
lution kernel, then combined into the reference 
frame by taking a median or average. After this 
was done, sources were detected and extracted on 
the reference frame using the SExtractor software 
( Bertin fc Arnoutslll996^ . When configured with 
a lower detection threshold, SExtractor extracts 
the number of spurious detections (e.g., diffrac¬ 
tion spikes around bright stars, or outer features 
of bright galaxies). These false detections were re¬ 
moved by careful visual inspection for each chip. 
The final catalog contains a total of 30,294 objects 


^ SWarp is a program that resamples and co-adds 
FITS images, distributed by Astrometic.net at 
http://www.astromatic.net/software/swarp. 


including both point and extended sources. 


3.1.2. Refining the Centroid of Each Object 


Prior to the photometry, the initial centroid co¬ 
ordinates of the target objects for each frame were 
computed b y using the WCSTools sky2xy routine 
( Mink 2002fl . The stored world coordinate system 
for each frame is used to convert the (R.A., decl.) 
coordinates from the master source catalog to the 
(x, y) pixel locations. However, actual positions 
of objects for each frame can be slightly moved 
from its original locations depending on the focus 
condition of instrument, seeing condition, and the 
signal-to-noise ratio (S/N) in the individual ob¬ 
servations. These types of positioning errors (i.e., 
centroid noise) will lead to the internal error for a 
photometric measurement that results from place¬ 
ment of the measuring aperture on the object be¬ 
ing measured. The situation gets worse for faint 
stars because the centroid position of them is itself 
subject to some uncertainty. The fractional error 
in the measured flux as a result of mis-centering is 
given by: 




1 A 2RA 

O’ 0-2 




( 1 ) 


where A is the positioning error, R is the radius of 
aperture, a is the profile width for a source with 
a Gaussian PSE, and total flux Fq (see Appendix 
A in Irwin et ah, 2007 for details). This expres¬ 
sion shows that if we set the aperture radius equal 
to the EWHM (i?=2.35-cr), even small differences 
in placement of the aperture (e.g., A=0.1-cr) may 
increase the uncertainty in the flux measurements 
(ft! 1 mmag). Thus, accurate centroid determi¬ 
nation is important to achieve the high-precision 
photometry. 

Eollowing the windowed centroid procedure in 
the SExtractor, a refined centroid of each object 
is calculated iteratively. On average, the rms un¬ 
certainty in the coordinate transformation using 
the WCS information was 0.044 ± 0.01 pixels for 
the bright reference stars. After the coordinate 
transformation from sky to xy, however, the cen¬ 
troid coordinates (xini, j/ini) are slightly misaligned 
from their actual ones (xcenter, 2/center)- The re¬ 
fined centroid values (xgnai, 2/finai) are used only if 
the maximum displacement is at least less than 1.5 
pixels. This condition prevents arbitrary shifting 
of a source centroid, especially for faint stars. 
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3.1.3. Estimation of Background Level 


1 


We estimated a local sky background by mea¬ 
suring the mode of the histogram of pixel values 
within a local annulus around each object, which is 
suitable choice for our uncrowded field (less than 
^1000 stars per chip). This process is a combi¬ 
nation of K-a clipping and mode estimation. The 
background histogram is clipped iteratively at ±3- 
(T around its median, and then the mode value is 
taken as: 

Mode = 3 X Median — 2 x Mean. (2) 

It represents the most probable sky value of a ran- 
domly chosen pixel in the sample of sky pixels 
( Stetsonl[l987ll . For relatively crowded regions, we 
utilized a background map created by SExtractor 
package using a mesh of 32 x 32 pixels and a me¬ 
dian filter box of 5 x 5 pixels. This map is used 
to conhrm the properness of individual sky values 
from annulus estimates. 


3.2. Multi-aperture Indexing Photometry 


Modern data reduction techniques aim to reach 
photon noise limit and minimize systematic ef¬ 
fects. For example, differential photometry tech¬ 
nique can be achieved b etter than 1% precision 


for brighter st a rs (e .g., lEverett fc Howelll 12001 


Hartman et al. 2005ll . and the deconvolution- 


based photometry algorithm leads to the min¬ 
imization of systematic effects in very crowde d 
fields (e.g., Magain et al.ll200^lGillon et al.ll2007 1. 
However, conventional data reduction methods of¬ 
ten fail to handle various artifacts in wide-held 
survey data. We present below a new photometric 
reduction method for precise time-series photom¬ 
etry of non-crowded helds, without the need to 
involve complicated and CPU intensive process 
(e.g., PSF htting or difference image analysis). 


3.2.1. Photometry with Multiple Apertures 

Our photometry is similar to standard aperture 
photometry, except in that we compute the hux in 
a sequence of several apertures and then determine 
the optimum aperture individually to each object 
at each epoch. This multi-aperture photometry is 
an efficient way to determine the optimum aper¬ 
ture size that gives the maximum S/N for a flux 
measurement. The maximum S/N is not neces¬ 
sarily at the same aperture for all objects, and 
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Fig. 2.— Tendency of optimal aperture selec¬ 
tion from the multi-aperture photometry. As stars 
get fainter, the optimal aperture decreases in size. 
Sources outside of this sequence are suspected to 
be contaminated. 


it ca n be obtained from a relatively small aper¬ 
ture ( Howell 1989ll . This photometric aperture is 
to achieve the optimal balance between flux loss 
and noises based on a relationship derived from 
the CCD equation (see Merline fc Howelll 1995). 
Figure [2] shows how the optimum apertures vary 
with the stellar magnitude. There is an obvious 
trend of decreasing aperture sizes with increasing 
magnitudes down to the faint magnitude limit in 
the example frame. 


Once we measure the flux of each object with 
the optimum aperture, we need to apply the aper¬ 
ture correction for small apertures. The aperture 
correction terms are estimated from the growth 
curve analysis of selected isolated bright stars (i.e., 
reference stars). The average curve-of-growth for 
each frame is calculated by measuring the differ¬ 
ence in magnitude between different pairs of aper¬ 
tures (up to 10 pixels aperture radius) and then 
an automatic correction is applied to all objects 
for each photometric aperture. The use of a com¬ 
mon aperture correction for each CCD assumes 
that there is no variation in the correction across 
the CCD. This flux correction method gives nearly 
the same brightness within the measurement un¬ 
certainties for all apertures. Any PSF variation 
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Fig. 3.— Example of multi-aperture photometry 
for one star (ID=10213 in the master source cata¬ 
log) through epochs tl to t4. Top panels: 100 x 100 
thumbnail images of the target star. Middle pan¬ 
els: normalized S/N as a function of aperture size. 
Bottom panels: aperture corrected magnitude as 
a function of aperture size. The arrows represent 
the peak locations in the aperture-S/N diagram 
(see text for details). 
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Fig. 4.— Typical example of multi-aperture in¬ 
dexing photometry. Top panel is r-hlter light 
curve of same star (ID=10213) as shown in Figure 
[3l Bottom panels show the multi-aperture index¬ 
ing scheme. The a;-axis is the aperture size and 
the y-axis is the differential magnitude between 
pairs of apertures Am(=miJ_ref — Wij.k)- We can 
see whether and at what aperture the differential 
magnitude (solid lines) begins to deviate from the 
model curve (dashed lines) for each epoch. 
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across the CCD causes systematic errors, however, 
and we deal with this in Section 3.4 and Section 

4. 


We performed the multi-aperture photometry 
based on the concentric apertur e photometry algo¬ 
rithm in DAOPHOT package ( Stetsonlll987ll . us¬ 
ing several circular apertures (up to 10 pixels aper¬ 
ture radius) with a fixed sky aunulus from 35 to 45 
pixels. The initial results of multi-aperture pho¬ 
tometry are stored in ascii-format photometry ta¬ 
bles, including the date of the observations (MJD), 
the pixel {x, y) coordinates, the aperture-corrected 
magnitudes with errors for each aperture, the sky 
values and its errors. Figure |3] shows the details 
of the multi-aperture photometry for one star at 
different epochs. In the former two epochs, the 
photometric apertures cau be properly selected by 
the S/N cuts, while in the latter two epochs, S/N 
increases for lager apertures. This unusual behav¬ 
ior is due to contaminatiou by a moving object. 
We automatically identifies similar unusual cases 
by the method of aperture indexiug. 


3.2.2. Determination of the Best Aperture with 
Indexing Method 

Our multi-aperture iudexiug method is simi¬ 
lar to the basic c oncept of the discrete curve- 
of-growth method ( StetsonI llQQOl l. Each object 
is indexed based on the difference in aperture- 
corrected magnitude between pairs of apertures 
Am(=mij_ref ~ ™ij,k) with mean trend for stars 
of similar brightness (see solid and dotted lines 
in bottom panels of Figure |4j respectively). The 
aperture with a 10 pixel radius is used as the fixed 
reference aperture. The mean trend is determined 
by computing the rms curve of the aperture cor¬ 
rection values for all measured apertures, and used 
to evaluate whether magnitude at a given aperture 
significantly differs from the mean trend. Since the 
rms value depends on the chosen magnitude inter¬ 
val, all stars are divided into groups according to 
their brightness in the individual frames. We de¬ 
termine the rms curve for each magnitude group 
using an iterative cr-clipping until convergence is 
reached. Objects lying within ±3 -ct of the model 
curve are indexed as contamination-free, and those 
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Fig. 5.— Histogram of differences in the magnitude offsets between the left- and right-side of each chip 
(Ai — Ar) for the whole data set. Each histogram is normalized by the total number of data frames N and 
is described by Gaussian distribution with different means and variances (dashed lines). 


above ±3-cr as a contaminated source. Figure |4] 
shows that multi-aperture indexing guides us to 
throw out some photometric measurements if they 
are discrepant from the mean trend. This ap¬ 
proach also gives us a chance to recover a measure¬ 
ment that would be otherwise thrown out. The 
problematic aperture can be simply replaced by 
one of the smaller apertures if it is indexed as 
contamination-free. This help us make a full use 
of the information offered by the data. 


3.3. Improved Photometric Calibration 

We present a new photometric calibration to 
convert the instrumental magnitudes onto the 
standard system, including a relative flux cor¬ 
rection of the left and right half-region of each 
CCD chip. As mentioned in the Section 3.1, 
MMT/Megacam shows the temporal variations 
in the gain between two amplifiers on each CCD, 
as well as between CCDs that are part of the same 
mosaic. It may have been caused by unstable bias 
voltage of the CCD output drain which has a 
profound impact on the gain of the output am- 
pliher. The level of readout noise is also unstable 
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between two amplifiers. To correct for this effect, 
the photometric calibration needs to be performed 
individually for each amplifier region. 

We use a sufficient number of (pre-selected) 
bright isolated stars as standard stars and com¬ 
pute the relative flux correction terms. These 
terms were derived for each frame using the mean 
magnitude offset (Aj j.) of standard stars {Ni^r) 
with respect to corresponding magnitudes (mjj) 
in the master frame chosen as an internal photo¬ 
metric reference 

Ai^r = — 'y ^ 

where m'ij is the aperture-corrected instrumen¬ 
tal magnitudes in other frames and ZP\ ,- is the 
photometric zero-points for the left- and right-side 
of each chip, respectively. To calculate the zero- 
points, we solve a linear calibration relation of the 
form: 


r - mij = - 0.07X -f 0.107(r - i), (3b) 

where r and i are standard m agnitudes from 


the ph otometric catalog of M3 7 (jHartman et al 
2008al) and X is an airmass term. The fit is per¬ 


formed iteratively using a sigma-clipping method. 

Figure [ 5 ] shows the difference in the magnitude 
offsets between the left- and right-side of each chip 
(Al — Aj.) for the whole data set, which is within 
±0.02 magnitude level for all 36 CCD chips. The 
histograms are normalized by the total number 
of data frames iV = 4, 730 and are described by 
a Gaussian function with slightly different mean 
values and shapes (dashed line). We clearly see a 
significant variation in difference between a pair of 
magnitude offsets for all CCDs. 


3.4. Field Distortion Correction 

The photometric calibration for wide-field 
imaging systems is also affected by position- 
dependent systematic errors d ue to a PSF vari¬ 


ation across the F OV (e.g., Ilvezic et al.l 12007 : 
Hodgkin et al. 2009l i. We derive PSF variations 


across the FOV with the SExtractor package. The 
change of the PSF shapes in the image plane is 
represented by spatial distribution of PSF FWHM 
values for several bright stars, with parameters of 
CLASS_STAR > 0.9, MAGERRJiUTO < 0.01 mag, and 
FLAGS = 0 (i.e., isolated point sources with no 



Fig. 6.— PSF FWHM variations as a function of 
distance from the image center for various seeing 
conditions. For each quadrant (denoted by differ¬ 
ent colors), the dashed lines represent the weighted 
spline approximation of the median value of each 
distance bin (one arcmin). 


contamination). Note that the PSF FWHM val¬ 
ues are defined as the diameter of the disk that 
contains half of the object flux based on a circular 
Gaussian kernel. Figure [6] presents the variation 
of the PSF FWHM as a function of distance from 
the image center for various seeing conditions. 
For each quadrant, the dashed lines represent the 
weighted spline approximation of the median value 
of each distance bin (1 arcmin). The result shows 
that the PSF FWHM varies significantly as a func¬ 
tion of position on the single-epoch image frames 
and variations are at the level of ^10% to 20% 
(0".l — 0".2) across the FOV. As the field distor¬ 
tion is not negligible from the center of field to its 
edges, such variations limit the accuracy of stellar 
photometry. 



















To address this issue, we perform a 2D polyno¬ 
mial fitting technique. For each frame, the correc¬ 


tion terms are described by a linear or quadratic 
polynomial depending on the position (x, y) only. 


J 


Amc.k(a;, y, merr)lin = Cok + ClkX + C2kfo 
A77ic,k(a;, y, merr)qud = Cok + ClkX + C2ky + C3kX^ + C4kj7^ + C5kX|/, 


where x, y are the pixel coordinates of N bright iso¬ 
lated stars, TTierr is the statistical weight in the fit¬ 
ting procedure, Ck are the sets of polynomial coeffi¬ 
cients for each aperture size, and ATOc,k(a;, y, nierr) 
are the difference in magnitude between the ref¬ 
erence aperture and k aperture, Amc^k(x, y) = 
hie, 20 ( 2 ;, 2/) — ^c,k(x,y), at the position (x,y) for 
each chip c. We derived the optimal parame¬ 
ter values from a nonlinear least-squares fit us¬ 
ing the Levenberg-Marquardt algorithm and au¬ 
tomatic differentiation!^ and choose between two 
models that best fit the data. 

Figure [7] shows the field-dependent magnitude 
offsets and the distortion correction by 2D polyno¬ 
mial fitting method for one example mosaic CCD. 
For the outer and the central region of the mo¬ 
saic, we compare the magnitude offsets between 
the reference aperture and the relatively smaller 
apertures as a function of (x, y) coordinates. Here 
x-axis is in the declination direction and y-axis 
is opposite to the right ascension direction. We 
find that the magnitude difference depends on po¬ 
sition (x, y) and is most discrepant in the outer 
part of the FOV. This effect is usually more sig¬ 
nificant in the y-direction than in the x-direction, 
especially for the case of aperture photometry per¬ 
formed with small apertures. The correction for 
field-dependent PSF variation reduces the initial 
^10% variation (gray lines) to less than ^1% 
(black lines). 

3.5. Photometric Performance Diagnos¬ 
tics 

3.5.1. Comparison with Archival Data 

We compare the photometric performance of 
the re-calibrated light curves with the non-de- 


®We used a LeastSquareFit module provided in the scien¬ 
tific python package (http;//www.scipy.org/ ). 


CCD 1 



r (mag) r (mag) 

Fig. 8.— Comparison of data recovery rate {top 
panels) and rms photometric precision of light 
curves {bottom panels) for the selected stars in the 
outer part (CCD 1) of the field of view. The black 
points are from the re-calibrated light curves, 
while the gray points are from the raw (left) and 
filtered (right) light curves in archive, respectively. 
Our results use nearly 100% of observed data and 
reach comparable accuracy without any filtering. 

trended archival light curvet@ by means of the two 
representative measures: (i) the rms photomet¬ 
ric precision Uph, and (ii) the data recovery rate 
-^recovery The former is defined as the standard 
deviation of light curves around the mean value as 


^Not e that the archi v al dat a have not been trend-filtered, 
but iHartman et al.l ll2009h used the de-trended archival 
light curves as part of the transit detection process in their 
paper IV of the prior analysis. 
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Fig. 7.— Example of the position-dependent magnitude offsets {gray symbols) and the distortion correction 
by 2D polynomial fitting method {black symbols) for selected CCD chips. For the outer {left) and the 
central {right) region of the mosaic, we compare the magnitude offsets between the reference aperture and 
the relatively smaller apertures (e.g., APS, APIO) as a function of {x, y) coordinates. The variation is usually 
more significant in the y-direction than in the ^-direction, especially for the case of aperture photometry 
performed with small apertures. 


a function of r magnitude: 


crph(iVj,mj) 


AT, 


\ P=1 


(wp - mj)2 
fVj-1 


(5) 


where Nj is the number of data points in each 
light curve, m-p is the observed magnitude, and 
rh-j is the mean magnitude of the object j, and 
the latter refers to the number of analyzed data 
frames normalized by the total number of observed 
data frames N for each object. In typical cases, 
the data recovery rate should be near unity in the 
bright magnitude regime and decreases with mag¬ 
nitude for fainter objects. For comparison, we de¬ 


cided to select light curve samples which show ei¬ 
ther no significant variability or seeing-correlated 
variations induced by image blending. We remove 
all known variable stars from the sample list based 
on a new catalog of variable stars in M37 field (S.- 
W. Chang et al. 2015, hereafter Paper II). To 
remove the light curves of blended objects, we use 
an empirical statistical technique to quantify the 
level of blending by looking for seeing-correlated 
shifts of the obje ct from its median magnitude 
( Irwin et al.ll2007l) . 

, _ X Xpoly 
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CCD 21 



r (mag) r (mag) 

Fig. 9.— Comparison of data recovery rate (top 
panels) and rms photometric precision of light 
curves {bottom panels) for the selected stars in 
the central part (CCD 21) of the field of view. 
The plotted symbols and notations are the same 
as Figure E] 


where 


X 


0-2 

P P 


( 6 ) 


for light curve points irip with uncertainties ap, 
and Xpoiy is the same statistic measured with re¬ 
spect to a fourth order polynomial in FWHM fit¬ 
ted to the data. We adopt the value b < 0.4 for 
the selecting light curves with no blending. The 
last selection criterion is that the light curves must 
exit both in the archive and our database. 

In the bottom panels of Figure |8] and Figure 
[HI we plot the rms photometric precision of light 
curves for the two Megacam CCD chips in the 
outer (CCD 1) and central (CCD 21) part of the 
FOV, respectively. The black points show the rms 
values of the re-calibrated light curves, while the 
gray points are for the raw and filtered light curves 
in archive. The first impression from this compari¬ 
son is that the typical rms scatter is overestimated 
from the raw light curves because of many out¬ 
liers in the photometric data (bottom left panels). 
For the better results, these light curves were fil¬ 
tered out in two steps: (i) clipping 5-cr outliers 


from each light curve and (ii) removing every data 
points that are outliers in a l arge number of light 
curves ( Hartman et ahilioOSbl l. In the second step, 
the outlier candidates are estimated by choosing a 
cutoff value for each CCD chip. The cutoff value 
is defined as the fraction of light curves for which 
a given image is a 3-cr outlier. This filtering was 
applied to remove bad measurements due to im¬ 
age artifacts or poor conditions, which were pre¬ 
viously thought to be unrecoverable, but resulting 
data loss is up to 20% of the total number of data 
points (bottom right panels). As shown in the top 
panels of the two figures, the data recovery rate 
for the re-calibrated light curves is close to 100% 
over a wide range of magnitude and it appears to 
be more complete compared with the raw (top left 
panels) and filtered (top right panels) light curves. 
At bright magnitudes (r < 16) the data recovery 
rate does not reach 100% because the exposure 
time was chosen to be saturated at a magnitude 
of r ~ 15. This comparison proves that our ap¬ 
proach is a powerful strategy for improving overall 
photometric accuracy without the need to throw 
out many outlier data points. 


Finally, we compare the light curves themselves 
for selected variable stars between the archive and 
our own. This comparison serves to illustrate how 
the photometric precision and data recovery rate 
of the time-series data affect the ability to address 
a variety of variability characteristics. Figure [101 
shows a direct comparison with the filtered light 
curves (top panels) and our re-calibrated light 
curves (bottom panels) for four variable stars. It is 
shown that our method recovers more data points 
(black) from the same data set of images. 


3.5.2. Comparison with PSF-fitting Photometry 


In order to further investigate possible system- 
atics in our approach, we conducted PSF-fitting 
photometry w ith DAOPHOT II and ALLSTAR 
(Stets^[lM3)- For each mosaic frame, we select 
bright, isolated, and unsaturated stars to make 
the PSF model varying quadratically with (x, y) 
coordinates. After PSF modeling, we run ALL- 
STAR to perform iterative PSF photometry of all 
detected sources in the frame with initial centroids 
set to the same values used for our own photom¬ 
etry. We then calculated apert ure corrections us¬ 
ing the package DAOGROW ( StetsoiJ 199fJl af¬ 
ter subtraction of all but PSF stars, which creates 
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Fig. 10.— Light curves for selected variable stars from Hartman et al. (2008b, top panels) and from this 
work (bottom panels), respectively. The recovered data points are marked with black squares. From upper 
left to lower right; ID=310010 (detached eclipsing binary star), ID=310139 (rotating variable star with 
flare), ID=100031 (pulsating variable star), and ID=100046 (aperiodic variable star). 


aperture growth curves for each frame and then 
integrates them out to infinity to obtain a total 
magnitude for each PSF star. The final step is 
to convert the instrumental magnitudes into the 
standard photometric system. For each frame, the 
initial zero-point correction is applied by correct¬ 
ing the magnitude offset with respect to the mas¬ 
ter frame. This places photometry for all frames 
on a common instrumental system. Following the 
same procedure in Section 3.3, the photometric 
calibration is performed individually for each am- 
pliher region. 

Figure [TT] shows the residual magnitudes and 
sky values between our multi-aperture photom¬ 
etry and the PSF-fitting photometry as a func¬ 


tion of position in the selected CCD chips. There 
are no position-dependent trends in the magni¬ 
tude residuals. For the brighter stars with r < 21 
mag, the rms magnitude difference between the 
two methods is very small (Achipi = —0.001±0.009 
and Achip 2 i = 0.001 ± 0.012, respectively), while 
for the relatively faint stars the rms difference is 
somewhat larger (Achipi = —0.002 ± 0.055 and 
Achip 2 i = 0.006 ± 0.050, respectively). The re¬ 
sults of this example indicate that we can reliably 
correct for the PSF variations by our calibration 
procedures. Meanwhile, our sky values are slightly 
higher than the ALLSTAR sky values, but not to 
the degree that can seriously affect photometric 
measurements. 
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Fig. 11.— Example of residual magnitudes and 
sky values between our multi-aperture photometry 
and the PSF-fitting photometry as a function of x 
and y coordinates, respectively, for selected CCD 
chips. The filled squares are relatively bright stars 
with r < 21 mag, while the open squares are stars 
with r > 21 mag. 


Figure [T^] shows a comparison of the rms dis¬ 
persion of the light curves obtained with our pho¬ 
tometry with respect to the that of the PSF-fitting 
photometry in the central region of the open clus¬ 
ter M3 7. We only compare the light curves of 
non-blended objects as described in Section 3.5. 
Our multi-aperture photometry does not reach the 
same level of precision as PSF-fitting photome¬ 
try for the faintest stars, while the PSF-fitting 
approach results in poorer photometry for bright 
stars. As shown in the right panels of Figure[T3 it 
is clear that our photometry tend to have smaller 
measurement errors with respect to the PSF pho¬ 
tometry for the bright stars. 



Fig. 12.— Comparison of the photometric preci¬ 
sion (rms) for our multi-aperture photometry and 
for the PSF-fitting photometry as a function of the 
r magnitude in the central region of the open clus¬ 
ter M37 (left panel). The arrows indicate the three 
different magnitude levels from bright to faint in 
our sample shown in the right panels. Note that 
the variable object IDs are taken from the new 
variable catalog of the M37 (see Paper II). 


(or unknown) systematics have been recognized 
by several exo-planet surveys because planet de¬ 
tection pe rformance can be easily damag ed b 


them (e.g. Kovacs et al.N2005l : iTamuz et al 


200 


I 


Pont et al. l2006[l . Also space-based time-series 


data (e.g., CoRoT and Kepler) are no excep¬ 
tion to this behavior although it is completely free 
from systematics caused by the turbulent atmo¬ 
sphere. Most of the raw light curves are affected 
by a secular (or a sudden) vari ation of flux without 


any obvious physical reason (iMazeh et al, 
iMislis et al.ll20 1 Otl Jenkins et al.ll20I0ll . 


20091: 


In order to check the properties of temporal sys- 
teinatics, we examined the correlation coefficients 
as measure of similarity between two light curves i 
and j obtained from a single CCD chip (CCD 1). 


4. TEMPORAL TRENDS IN THE RE¬ 
CALIBRATED LIGHTCURVES 


. _ 1 (^) nLjLj 

'j “ A - 1 CTiCTi ’ ^ ’ 


From a visual inspection of the re-calibrated 
light curves in the same CCD chip, we found 
that some light curves tend to have the same 
pattern of variations over the observation span 
(Figure IT^ . This kind of systematic variation 
(i.e., trend) is often noticed in other studies. For 
example, the importance of minimizing known 


where L{t) is the flux of each star at time t, N is 
the total number of measurements, L is the mean 
flux of each star, and a is the standard deviation of 
L{t). This comparison is a point-by-point compar¬ 
ison and is done for every pair of light curves in the 
data set. The resultant similarity matrix can be 
used to identify correlated pairs of light curves and 
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Fig. 13.— Example of shared systematic trend 
in the light curves of isolated bright stars. The 
numbers on the upper left- and right-side of each 
panel show the star’s identification number and 
the mean value of magnitude, respectively. The 
numbers on the s-axis indicate the corresponding 
timestamps in each frame. The deviations from 
the mean value, Ar, are less than ±0.01 mag level 
(with a rms values of 0.0021, 0.0025, and 0.0033 
mag from the top panel down). These stars show 
a similar pattern of light variations over the ob¬ 
servation span. 


to determine whic h light curve is least like all other 


light curves (e.g., Protopanas et al. 200dl . After 


that, we selected stars showing most systematics 
based on a hierarchical cluste ring method with the 


et al. 20091 for 


correlation coefficients (See 
more details). Figure [14] represents spatial dis¬ 
tribution of the most prominent trend groups on 
the CCD plane (top panel) and its strongly cor¬ 
related features determined by the weighted sum 
of normalized light curves (bottom panel). There 
are two interesting features in this figure: the hrst 
one is that each trend covers only a certain part 
of the sky area and the second one is that some 
portions of neighboring trends show different vari¬ 
ation patterns even at the same moment in time 
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Fig. 14.— Top panel: spatial distribution of 
the most prominent trend groups on the CCD 
plane which covers the small areas of the sky 
(6'. 14 X 2'.73). The colored dots represent the 
most systematic stars classihed as group 1 (red), 
group 2 (blue), and group 3 (green), respectively. 
One interesting property of trends is that they are 
localized within a CCD frame where stars are iso¬ 
lated from each other. Bottom panel: systematic 
features extracted from the selected light curves. 
These are calculated by the weighted sum of nor¬ 
malized light curves from each group. 


(shaded gray region in the hgure). In particu¬ 
lar, we found an anti-correlated variation for the 
trends between the group 1 (gi) and the group 
2 ( 32 ), so we might expect to hnd possible noise 
sources that are responsible for these discrepan¬ 
cies. Why the trends are different and localized 
within a single CCD frame is a subject of further 
study, but it is probably related to subtle changes 
in point spread function and sky condition within 
the detector FOV. 

For these two groups, we consider a possible 
causal relationship between the systematic trends 
and average object/image properties. Figure [T51 
shows the differences in trend, differential mag¬ 
nitudes, sky level, and PSF FWHM between the 
two groups, respectively. In the top panel, we plot 
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Fig. 15.— Example subset of data for comparing 
the trends with several parameters from the two 
groups {gi and 52 ) of stars. The numbers on the 
cc-axis indicate the corresponding timestamps be¬ 
tween the frame 1 and the frame 400. From the 
top to bottom, the panels show the average dif¬ 
ference of trend, flux concentrations in which we 
applied a —0.02 mag shift to the m 2 o — mi 2 value, 
sky level, and PSF FWHM value between the two 
groups. 


the magnitude difference in trends, which shows 
variations in the range of ±0.02 mag. We suspect 
that this may be due to the different concentra¬ 
tion of star light between these two groups. It 
can be checked by using the magnitude difference 
Am = m 2 o — mop, where m 2 o and map are the 
reference aperture and the relatively small aper¬ 
ture, respectively. In fact, we already know that 
there is a magnitude variation in Am depend on 
the aperture size due to the field-dependent PSF 
variation (see Section 3.4). For example. Figure 
1 161 shows the response of multi-aperture photome¬ 
try for the two groups of stars at one epoch (MJD 
= 53726.14817) before and after applying the dis¬ 
tortion corrections. Although the magnitude vari¬ 
ation between the group 1 and the group 2 seems 



Fig. 16.— Systematic variations in the magnitude 
offsets. Am, for the original multi-aperture magni¬ 
tude measurements as a function of aperture size, 
which is marked with a gray dots in all panels. A 
positive Am indicates that the photometric mea¬ 
surements with the corresponding aperture size, 
map, are brighter thau those of reference aperture 
m 2 o, while a negative Am indicates vise versa. 
The dashed lines are the rms model profiles in¬ 
troduced in Section 3.2.2. There is a noticeable 
distinction between the group 1 {top panels) and 
group 2 {bottom panels) when looked at differ¬ 
ent concentration levels. But the treuds are not 
fully explained by the difference patterns in Am 
because those correlated variations become small 
after the correction for the PSF variation (black 
dots). 


to have different behavior as a function of aperture 
size (gray points), it is negligible after the removal 
of the PSF variation (black points). 

We also check that the possible contribution 
of sky level (Asky) and PSF FWHM differ¬ 
ences (AFWHM) to the systematic trends on 
the re-calibrated ligh t curv es. As mentioned by 
Bramich fc Freudlind ( 20I2[l , sky over-subtraction 


may lead to the systematic trends as a function of 
the PSF FWHM, the amplitude of which increase 
for fainter stars. The third and forth panels of Fig¬ 
ure show the variation of the mean Asky and 
AFWHM, respectively. In our case, however, the 
form and amplitude of trends seem independent 
of sky level and PSF FWHM. 

Some other possible sources that may con¬ 
tribute to the observed systematic trends include: 
higher order variations in the PSF shape beyond 
just the FWHM; cross-talk from other amplifiers. 
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or ghosts from bright stars undergoing multiple re¬ 
flections within the optics; non-uniform variations 
in the gain; and unmodeled temporal atmospheric 
variations that are dominated by Rayleigh scatter¬ 
ing, molecular absorption by 
por, and aerosol scattering 
2008h . While we find a clear presence of trends 
that should be removed, we are not able to iden¬ 
tify their exact cause. 


ozone and water va- 
(Padmanabhan et al. 


4.1. Removal of Temporal Systematics by 
Photometric De-Trending (PDT) 

In order to reduce systematic effects in photo¬ 
metric time-series d ata, several method s were in- 
troduced (e.g., TFA : Kovacs et al ]|2005t Sv s-Rem: 


^muz_e^^L 20^; PDT: Ki^_e£_^ J200^_CDA: 


Mislis et al.l 201(1 PPG: Twicken et al. 20101) . All 


of these algorithms share a common advantage 
that they work without any prior knowledge of the 
systematic effects. We use the PDT algorithm, 
which has been designed to detect and remove 
spatially localized patterns. By default, this al¬ 
gorithm works with a set of light curves that con¬ 
tain the same number of data points distributed in 
the same series of epoch. In many cases, however, 
missing data occur when no photometric measure¬ 
ments are available for some stars in a given ob¬ 
served frame. These missing data can be simply 
replaced by means, medians, or the values from the 
interpolatio n of adjacent data points in each ligh t 


interpolatio n ot ad.iacent data points in eacn iign t 
curve (e.g., iKovacs et al.ll2005t iKim et all 120101) . 


Although using the replaced value is the easiest 
way to reconstruct the light curve to be analyzed, 
it is not appropriate if the time separation between 
two subsequent observations is too large. Instead 
we use more straightforward approach by apply¬ 
ing the PDT algorithm in two separate steps: (i) 
we construct the master trends from the subset 
of bright stars, and (ii) de-trend light curves of all 
stars with most similar master trend and matching 
time line. 

We briefly describe the main procedure of our 
de-trendi ng process follow ing the algorithm de¬ 
rived by iKim et all (j2009[ ). We first select the 
template light curves from bright stars that show 
the highest correlation in the light-curve features. 
The total length of template light curves should 
be long enough to cover the whole time span of 
observations. In this step, we take a sequence of 
data points, as the reference time line. Us- 
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Fig. 17.— Example of selected light curves be¬ 
fore (gray) and after (black) the removal of sys¬ 
tematic trends. The numbers on the x-axis in¬ 
dicate the corresponding timestamps between the 
frame 1 and the frame 2000. The y-axis is r-filter 
magnitude (normalized by its mean value). While 
the morphology of the two light curves in the left 
panels appear to be variable stars of some kind, 
these turn out to be non-variable after applying 
the photometric de-trending method. In the case 
of the right panels, all true variabilities are pre¬ 
served from the raw light curves. From upper left 
to lower right: ID=10032, ID=10039, ID=170088, 
and ID=170108. 


ing the correlation matrix calculated from Equa¬ 
tion (7), we extract all subset of light curves that 
show spatio-temporally correlated features (i.e., 
clusters). Each cluster is determined by hierar¬ 
chical tree clustering algorithm based on the de¬ 
gree of similarity. Next, we obtain master trends 
Tc(tref) for each cluster by weighted average of the 
normalized differential light curves, /i(tref): 


T'c(Uef) = 


EiJl WiMtref) 

Ei=i 


/*(Uef) = 


Aj(tref) 7/j 


t2 ’ 
^fi 


( 8 ) 
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where Nc is the total number of light curves in 
each cluster c, Li is mean value of ith light curve, 
and af- is the standard deviation of /i(tref)- After 
determining the master trends, we de-trend the 
light curves of all stars with matching master trend 
and time line. We adjust the temporal sequence of 
measurements for the master trends Tc(ti) by that 
of individual light curves to be de-trended Li{ti). 
Because each light curve is assumed as a linear 
combination of master trends and noise, we can 
determine the optimal solution by minimizing the 
residual between the master trends and the light 
curve: 


h{U) = '^l3icTc{ti) +ei{ti), 

C=1 


MU) = (9) 

where m is the total number of master trends and 
Pic are free parameters to be determined by means 
of minimization of noise term eptM. 

Figure [T7] shows examples of our light curves 
before and after removing the systematic trends. 
The algorithm we used for de-trending removes 
only the systematic variations that are shared by 
light curves of stars in the adjacent sky regions 
(left panels), while all kinds of true variabilities 
are preserved (right panels). 


5. IMPACT OF THE NEW CALIBRA¬ 
TION OF PERIOD SEARCH 


The usefulness of our photometry is tested for 
a set of variable stars. We immediately find abun¬ 
dant cases of improvements in the following three 
aspects: (i) refinement of the derived period, (ii) 
detection of a new significant peak in the pe- 
riodogram, and (iii) separation of non-variable 
candidates where systematics in the light curves 
were mistaken for true variability. For each case, 
we compare light curves and power spectra for 
archival data and our data. 


For the first example, we show that a new pho¬ 
tometric measurement and calibration allowed us 
to derive a much improved refinement of the light 
curves and of the derived periods iFigure ITSl) . 
We performed a Lomb-Scargle (L-S: Scarglelll982l) 
search of both archival and new light curves for 
periodic variable star (V427). The light curves 
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Fig. 18.— Top panels: archival (left) and our fi¬ 
nal (right) light curves of a periodic variable star, 
V427. The light cur ves are folded by the best- 


fit period of 5.4615 ( Hartman et al. 2008bll and 


4.4158 days, respectively. Such an period differ¬ 
ence for the same star comes from data itself (e.g., 
both the different noise levels and the different 
data sampling intervals). Bottom panels: the re¬ 
sulting amplitude spectrum of each light curve is 
calculated with PERI0D04 package. 


are folded by the best-fit period of 5.461(0 and 
4.4158 days, respectively. We also calculated the 
false-alarm probability (log FAP) for each peak 
and its signal to noise ratio (S/N): log FAP archival 
= —28.37, S/Narchivai = 43.2 for the archival data 
and log FAPnew = -181.07, S/N„ew = 84.5 for 
the new data. Since we can get better estimation 
with much lowered minimum FAP value, our new 
period is the most likely result. In the bottom 
panels, the resulting amplitude spectrum was cal¬ 
culated with PERI0D04 package. Since the archival 
data is more noisy than the new one, it is rather 
complicated to interpret the peaks of its power 
spectrum. 

For the second example, we show the newly 
discovered low-amplitude pulsating variable star 
(Figure [121). We used the PERI0D04 package to 
find multiple pulsation periods. The whole pro¬ 
cess of identifying, fitting, and pre-whitening suc¬ 
cessive frequencies was repeated until no signifi- 


^For the archival data, w e adopt the period fou nd on the 
filtered light curves from iHartman et al.l ll2008bf ). 
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Fig. 19.— Top panels: archival (left) and our fi¬ 
nal (right) light curves of low-amplitude pulsating 
variable star, V2276. The archival data is not ade¬ 
quate to discriminate a signal of astrophysical ori¬ 
gin from the noise of the data stream. Bottom pan¬ 
els: the resulting amplitude spectrum of each light 
curve is calculated with PERI0D04 package. A sig¬ 
nificant peak is detected at the frequency 22.3979 
days“^ (indicated by arrow). The phased diagram 
of this candidate frequency is shown in the low- 
right panel. We see that the model fits well the 
overall pulsating variability (red line). 


cant frequencies were found. We adopt a conser¬ 
vative approach in selecting the statistical signifi¬ 
cant peaks from the amplitude spectrum. A S/N 
amplitude ratio of 4.0 is a good criterion for inde¬ 
pendent frequ encies, equivalent t o 99.9% certainty 
of variability ( Breeer et al.lll9f^ . While no clear 
periodicity was found in the archive data, our am¬ 
plitude spectrum shows a clear excess of power 
centered at 22.3979 days“^ with peak amplitudes 
of about 1 mmag (S/Nnew = 9.02). 


The last example is the opposite case of the sec¬ 
ond. Figure [201 shows that this object is unlikely to 
be a variable source because there is no evidence 
for any significant peaks, which indicates that the 
variations are mostly noise. Extensive study on 
variabilities will be presented in Paper II. 


6. CONCLUSION 

In this paper, we introduce a new time- 
series photometry with multi-aperture indexing 


Fig. 20.— Top panels: archival (left) and our final 
(right) light curves of a star, V347. This star turns 
out to be non-variable in the new data. Bottom 
panels: all other details are same as Figure [TOl 
but the subfigure of low-right panel is 200x200 
thumbnail image of the target star. There is no 
potential sources of contamination to hamper the 
interpretation of the power spectrum. 


and spatio-temporal de-trending techniques, to¬ 
gether with complex corrections to minimize in¬ 
strumental biases. We used the archival, high- 
temporal time-series data from one-month long 
MMT/Megacam transit survey program. The re¬ 
calibration of the archival data has made several 
improvements as follows: (i) the photometric in¬ 
formation derived from the multi-aperture index¬ 
ing measurements is useful to obtain the best S/N 
measurement, but also to diagnose whether or not 
the targets are contaminated; (ii) the resulting 
light curves utilize nearly 100% of available data 
and reach precisions down to sub mmag level at 
the bright magnitude end without the need to 
throw out many outlier data points, which makes 
it possible to preserve data points that show in¬ 
trinsic sudden variations such as flare events; (iii) 
corrections for position-dependent PSF variations 
and de-trending of spatio-temporal systematic 
trends improve the quality of light curves; and 
(iv) new photometry enables us to determine the 
variability nature and period estimate more accu¬ 
rately. 

While this study deals with a particular set 
of data from MMT, we find our approach has a 
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potential for other wide-held time-series observa¬ 
tions. Multi-aperture indexing measurement is a 
powerful tool in isolating and even correcting var¬ 
ious contaminations. Spatio-temporal de-trending 
is also very useful in removing systematics caused 
by PSF variation and even non- uniform extinc¬ 
tion o f thin clouds across the FOV. IChang fc Bvun 
( 20 Fil l proved this for different sets of archival sur¬ 
vey data. 
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